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Abstract 

We study the decay processes B — * Kir by using generalized decay amplitudes with both the final 
state re-scattering strong phases and the topological quark diagrammatic strong phases included 
together as part of fitting parameters. Using a generalized approach, so called "isospin phased 
topological approach", for all the currently available data of B — > Kir decays, we determine the 
allowed values of the relevant theoretical parameters, corresponding to the electroweak penguin, 
the color-suppressed tree contribution, strong phase differences, etc. In order to find the most likely 
values of the parameters in a statistically reliable way, we use the x 2 minimization technique. We 
find that the long distance final state re-scattering, when taken at proper value, can provide a 
reasonable fit to the standard model with the perturbative QCD estimated values, and therefore, 
it is premature to conclude that it requires new physics to explain the CP violating B — > Kir data. 
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I. INTRODUCTION 



There are four different decay channels (and their anti-particle decay channels) for B — > 
Kit processes, depending on the electric charge configuration: B + — > K°tt + , B + — > K + n Q , 
5° — > K + tt~, and 5° — > K°tt°. All the S — > i^7r modes have already been observed 
in experiments and their CP-averaged branching ratios have been measured within a few 
percent errors by the BaBar and Belle collaborations . The observations of 

the direct CP asymmetry in B° — > K ic* have also been recently achieved at the 5.7 a level 
by BaBar and Belle [zl la S Q> C I • For the other B — > Kir modes, the experimental results 
of the direct CP asymmetries still include large errors. Certain experimental data (e.g., 
the branching ratios (BRs)) for B — > Kit are currently more precise than the theoretical 
model predictions based on QCD factorization (QCDF), perturbative QCD (pQCD), and so 
on. Thus, these decay modes can provide very useful information for improving the model 
calculations, and at the same time, the model-independent study becomes very important. 

In the light of those new data, including the direct CP asymmetry in B° — > K ± ti t , 
many works have been done to study the implications of the data 

Q Q, 20, 21 1 . The quark level subprocesses for B — > Kit decays are b — > sqq (q = u,d) 



penguin processes, which are potentially very sensitive to any new physics effects beyond the 
standard model (SM). Thus, with the currently available precision data, it is very important 
to investigate these modes as generally and critically as possible. 

At the bottom of this important problem, lies the fact that final state re-scattering phases 
by strong interaction play a crucial role in generating CP violation. Part of the strong 
interaction phases (as short- distance strong phases (sdSPs) related to topological quark 
diagrams) can be investigated by using various QCD models such as QCDF or pQCD, even 
though the results are quite model dependent. However, there are also long-distance final 
state strong phases (IdSPs), that are very difficult to calculate due to hadronic interactions at 
low energy scale (even thought some attempts had been done). Despite facing the fact that 
we do not understand strong phases very well and that there are already many theoretical 
claims of possible new physics from the data of B — > Kir, here we are trying to be as model- 
independent as possible in fitting the B — * Kit data, by choosing both sdSPs and IdSPs as 
free fitting parameters, to see if new physics is still required by the recent data or not. 

In this work, we study the decay processes B — > Kir by using generalized decay ampli- 
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tudes within the so-called "isospin-phased topological approach", i.e. with both the strong 
phases of the final state re-scattering (IdSPs) and the strong phases of the topological quark 
diagrammatic origin (sdSPs). If we ignore topological strong phase differences, this general 
approach becomes the original isospin approach. Inversely, if we ignore the strong phase 



221, which is 



differences in isospin amplitudes, it reduces to the normal top olog ical analysis 
based on an exact flavor SU(3) symmetry in B meson decays 231 ] . 

We note that the importance of final state interactions (FSI) has recently been recognized 



again in hadronic B decays 



24j . in which the importance of the existence of the soft final 



state re-scattering effects has been pointed out, esDecially in a model-independent topolog- 
ical quark diagrammatic analysis of B — > Dn |25| and B — > Kit decays. It is certainly 
conceivable that FSI can modify the predictions based on the short distance diagrammatic 
analysis. Likewise, the branching ratios of certain decay modes B s ^ — > are expected to 
be extremely small if re-scattering effects through FSI do not alter the predictions of the 
diagrammatic approach. We notice that the soft final state re-scattering effects can violate 
the exact flavor SU(3) symmetry even though the isospin symmetry will still hold with FSI, 
since FSI may rather be at low energy scale of light final particles. In other words, at the 
scale that FSI are activated, the flavor SU(3) is broken but the SU(2) isospin symmetry is 
still valid. And FSI can be parameterized as the isospin phases in the limit of the elastic 
re-scattering. 

We note that including both sdSPs and IdSPs may amounts to a possible double counting: 
Since they typically involve physics at different scales with different symmetries like flavor 
SU(3) symmetry or SU(2) isospin symmetry, the double counting may not be as serious as 
one thinks. As is well known, the scattering of hadrons exhibits a two-component structure 
of "soft" and "hard". We associate the high scale hard scattering components between 
pointlike constituents with the SU(3) strong phases of quark diagramatic origin, and the 
low scale soft components of FSI with the SU(2) isospin phases. In addition, since they are 
just taken as model independent fitting parameters, even if there are some double counting 
of physical effects contained in our fitting parameters, it is not going to affect our final 
physical conclusions regarding new physics. 

Here we are mainly interested in investigating whether the conventional SM predictions 
are consistent with the current data even after we include FSI effects. Furthermore, if 
there are some deviations between the conventional estimates and the experimental results, 
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we intend to identify carefully the source of the deviations and estimate how large the 
contribution from the source can be. Then, by comparing our result with the conventional 
SM predictions, we shall be able to verify whether the current data indicate any new physics 
effects. In order to find the most likely values of the theoretical parameters in a statistically 
reliable way, we will adopt the x 2 analysis. 

We organize the paper as follows: In Section II, we introduce all the relevant formulas 
for B — > Kir decays, step-by-step generalizing the decay amplitudes. We also give formulas 
for BRs, direct CP violations and mixing-induced direct CP violation. In Section III, we 
present the summary of the recent experimental results on B — > Kir modes. And we do % 2 
analysis for all the experimental observables, and discuss its physical implication within our 
generalized approach. In Section IV, we give conclusions. 

II. THE DETAILED FORMULAS FOR B -> Kir DECAY MODES IN THE 
ISOSPIN-PHASED TOPOLOGICAL APPROACH 

We first introduce the decay amplitudes within the generalized "isospin-phased topolog- 
ical approach", then we summarize the formulas for the relevant decay amplitudes, BRs, 
direct and indirect (mixing-induced) CP asymmetries of B — > Kir processes. 

A. Fusion of Isospin and Topological Approaches 

We write 73-meson decay amplitude, e.g. B° — > ir~K + , by including both the topological 
strong phases and the isospin related strong phases, so called "isospin-phased topological 
approach", as: 

a = J2 A i e M^i) (i) 
i 

= J2i D i exp(i(Ji) + D 2 exp(i8 2 ) + exp(iA 7 ), (2) 
i 

where Aj is an isospin amplitude with the isospin related (IdSPs) strong phase Aj, and we 
decompose each isospin amplitude into a sum of the topological diagrammatic amplitudes, 
Di, with the topological strong phases (sdSPs) Si. As can be seen below, for B° — > ir~ K + 
decay channel Aj's are expressed with 3 isospin amplitudes A 3 / 2 , Ay 2 and B ± / 2 , with strong 
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phases a' 3 / 2 , ct[/ 2 and (3[/ 2 , respectively. And each Aj is expressed with a sum of the topolog- 
ical amplitudes, T, C, P, Pew etc. with strong phases 5' T , 5' c , 5' P , 5' EW etc., respectively. Now 
if we ignore phase differences in the topological amplitudes, we recover the original isospin 
amplitude with a positive real value \Aj\ (except for the relevant weak phases). Inversely, if 
we ignore phase differences among isospin amplitudes, the framework reduces to the normal 
topological analysis. 

As final states, we take ir + = ud,ir° = dd ^ u , 7r~ = —du, for pions and K + = us,K° = 
ds, K° = sd, K~ = —su, for kaons, and for B mesons B> + = bu, B° = bd, B° = bd, B~ = —bu. 
Isospin relations are 



A(B° - 




= X = 


A3/2 + A1/2 — B1/2, 


(3) 


V2A(B° - 


- n°K°) 


= Y = 


2A3/2 — A1/2 + B1/2, 


(4) 


A(B + - 




= Z = 


A3/2 + Ay 2 + B1/2, 


(5) 


V2A(B + - 


- n°K + ) 


= w = 


2^3/2 — A1/2 — B1/2, 


(6) 



where A and B represent AI = 1 and AI = components, respectively, and subscripts 
denote the isospin of final states. These quantities satisfy the isospin relation, 

X + Y = Z + W. (7) 

Decay amplitudes can be written in terms of topological contributions 

A(B° -> n-K + ) <-> -P - P% w - T, (8) 

A(B° - *°K°) ~ ±P - j= 2 P EW - j= 2 C, (9) 

A{B + ^7i + K°) <-> P + A, (10) 

A(B+ - J>K+) - -j-P - ±P EW - j= 2 P£ w --feT-j-C- j= 2 A, (11) 

where T, P, Pew-, C, P E wi ^Ci A are (T) tree, (P) penguin, (Pew) electroweak (EW) pen- 
guin, (C) color suppressed tree, (P E w) color-suppressed EW penguin, and (A) annihilation 
amplitude, with strong phases S' T , 5 P , S' EW , 5' c , etc., respectively. P represents the combina- 
tion P = P'~ \Pgy. 

Now we introduce isospin phase 0:3/2 > a 'i/2 an d fl'1/2 ^° eacn isospin component 
(A 3/ 2,A 1/ 2,B 1/2 ), we get: 

A 3/ 2 = \{~T -Pew-C- Pgy)e^, (12) 



A 1/2 = \{-T + 2P EW + 2C- P% w + 3A)e ia 'v*, (13) 
o 

Bi/2 = \(2P + T + P% W + A)e l/3 '^, (14) 

where we ignored A (weak annihilation contribution) because it is estimated to be much 
smaller than others within the SM evaluations of pQCD and QCD factorization. 



B. Branching Ratios, Direct and Mixing-Induced CP Asymmetries 

In this subsection we express branching ratios, direct and indirect CP asymmetries. We 
remark that there exists a (naive) conventional hierarchy within the SM among the topo- 
logical diagrammatic contributions: 

1 > r T ~ T EW > r C ~ T EW > T 'a > ( 15 ) 

where 

r T = \T\/\P\, r BW = \P EW \/\P\, r a = \C\/\P\, = \P$ W \/\P\, r A = \A\/\P\. 
For instance, in the pQCD approach, those ratios are roughly estimated as 12 



r T w 0.21, r EW w 0.14, r c w 0.02, r C Ew w 0.01, r A w 0.005 . (16) 

It is also known that within the SM under flavor SU(3) symmetry, the relation 5 T ~ 5 EW 
holds to a good approximation [27( , which can be deduced from the fact that the topology of 
the color-allowed tree diagram is similar to that of the EW penguin diagram. Here we neglect 
the tiny quantities r A and . However, because recent studies on two-body hadronic B 
decays show that the color-suppressed tree contribution could be enhanced to a large amount 
through certain mechanisms j3, 13 13 we ^ ee P in order to take that possibility into 



12, 291, where all the linear terms for 



account. This treatment differs from that in Refs. 
r A and r^ w as well as r c were simply neglected. We note that the physical strong phases, 
which appear in the branching ratio and CP asymmetries after taking out overall phases, 
are defined to be 

5t = 5't — 5'p, ^ew = $'ew ~ ftpi &c = ft c — ft P) 5 EW = ft EW — ftp, 

and 



"1/2 = a'i/2 - P[/2 and "3/2 = "3/2 _ • 
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T T T, 



a. Branching ratios The CP-averaged branching ratios are given by, 

B + - = B(B d -> KV) 

oc |P| 2 |A,| 2 {l + 2r T cos0 3 ^ -2r w Cf w 
+r 2 T \Ax\ 2 + r 2 EW \X x \ 2 - 2r T r EW cos(0 3 )C* 
2B +0 = 2B{B ± -> ^tt ) 

oc |P| 2 |A t | 2 {l + 2r T cos0 3 ^-2r BM/ ^ 

+r 2 |A W | 2 + r 2 JE W | 2 - 2r T r EW cos^OC^ + 2r c C^} , 
£ 0+ = S( J B ± -> Ktt*) 

oc |P| 2 M 2 {l -2r T cos0 3 C r Z T +2r w Cf w 
+r 2 \A Z \ 2 + r 2 j£ z | 2 - 2r T r EW cos(0 3 )C£ 
2B 00 = 2B[B d -> Ktt ) 

oc |P| 2 |A t | 2 {l - 2r T cos0 3 C r Y T + 2r w C r Y w 
+rJ|A Y | 2 + r 2 J£ y | 2 - 2r T r w cos(0 3 )^ 

Here for simplicity we have defined parameters listed below: 



-iY 



T 

' T 

X 



C. 



= Re 
= Re 
= Re 
= Re 



e i&T A x 
ie i5c E x 
ie i&T A x 

e i{5T - 5Ew) A x ^x 



C 



'X — 



, S. 



' c 

x 



= Re 

= Re 
= Re 



e iS E W£ x 



ie i8c H x 
ie l(5T ~ 5Ew) A x i: x 



where 



Ax 


= e 


-^1/2 


V 2 


Ay 


= e 


-^1/2 




A z 


= e 


"^ /2 




A w 


= e 


"^/2 





2 



6 

1 .V 





e -^ /2 


Sy = 


e -^/2 






S w = 


Sy. 



1 

3 f 



(17) 



(18) 



(19) 



(20) 



2 



1 ictf, ,„ 1 ia' 



3 6 

The CKM elements are A t = V t * b V ts , \ u = V£ & V^ S and </> 3 (= 7) is the angle of the unitarity 
triangle. 

b. Direct CP asymmetries The direct CP asymmetries are given by 



+ _ , B(B° -> K-7I+) - B{B° -> ir+71-) _ |X| 2 - |X| 
CP ~ £(P° -> K-7T+) + B(B° -> K+7T-) ~ |X| 2 + |X| 
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A+0 
J\ C P 



-2r T sm <f> 3 S* + 2r' T sm(2<f) 3 )S*C* 

+2r T r EW sin hSr T r EW ~ ^ r T r E w sin ^ S r T C r EW - 2r c sin hSf c , (21) 
B(B~ -> ^-tt°) - ff(ff+ -> if +tt°) _ |T¥| 2 - |V^1 2 
-6(5- -»• ^-7T°) + £(5+ -»• ir+vr°) ~ |W/| 2 + \W\ 2 
-2r T sin + 2r 2 sin(20 3 )^Cf 



+2r T r w sin faS™. - ±r T r EW sin 3 S£ - 2r c sin 3 S r ^ , (22) 



B(B~ -> ^°7T-) - £(£+ -> K°7r + ) 




2 _ 


Z| 


2 


<B(B- -> K°7T-) + £(£+ -> K°7T+) 




2 + 


Zj 


2 



+2r T sin <I) 3 S? t + 2r 2 T sin{2^)S^C^ 

+2r T r EW sin hSr T r EW ~ 4r T r w sin fo S r T Cr EW + 2r c sin <p3 S r c , (23) 
B(B° ^ K\°) - B(B° ^ K°ir°) \Y\ 2 - \Y\ 2 



00 

CP ~ B(B° -> i?°vr ) + £(5° -> iW) |?| 2 + |y| 2 
~ +2r T sin03<+2r2sin(20 3 )<C r y T 

+2r T r EW sin 03^ w - 4r T r EW sin 3 ^ C^ w + 2r c sin 3 S^ . (24) 

Notice that considering the conventional hierarchy given in (|15p and (|16|1. the direct 
CP asymmetries Ac°p (j22j) and -4jp (}2*T]) are expected to be almost the same including 
their signs, because the dominant contribution to them is identical. However, the current 
experimental data show that A^° P and A^p are quite different from each other and even 
have mutually opposite signs, as shown in Table HJ 

c. Mixing-induced CP asymmetry The time-dependent CP asymmetry for B° — > K s tt° 
is defined as 

T(B°(t)^Ky)-T(B°(t)^Ky) 



A Kc; no(t) = 



T{B°{t) -> K S 7T°) + T{B°{t) -> K S 7T°) 

S K o sin(Am d t) + C K o cos(Am d t), (25) 



where T denotes the relevant decay rate and Am^ is the mass difference between the two B° 
mass eigenstates. S Ks7T o and C K w o are CP violating parameters. In the case that the tree 
contributions are neg lected for B° -> K s ir° } the mixing-induced CP violating parameter 
Sk s tt° is equal to sin(20i) [0i (= (3) is the angle of the unitarity triangle]. The expression 
for Sk s tt0 (up to r order) is given by 

~ sin 20! - 2r T sin 3 Re ( e - 2 *e*A y ) - 2r c sin 3 Re (e^^'e^Sy) . (26) 
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The measured value of Sk^o (Table |l} is different from the well-established value of 
sin(20i) = 0.725 ± 0.037 measured through B -»• J/tpK^ [30]. It may indicate that the 
subleading terms including r c and r EW in Eq. (J26|) play an important role. In Ref. j^lj, the 
authors showed that as Sk s k° varies, the allowed region for r c varies very sensitively, but 
that for r w not, as can be seen from Eq. (|26|). 

III. THE B -> Kir PUZZLE AND ITS PHYSICAL IMPLICATIONS 
A. Summary of Present Experimental Results 

We first summarize the present status of the experimental results on B — > Kir modes 
in Table I, which includes the BRs, the direct CP asymmetries (Acp), and the mixing- 
induced CP asymmetry (Sr^o). We see that the averages of the current experimental values 
for the BRs include only a few percent errors. Furthermore, the direct CP asymmetry in 
B° — > K ± -k t has been recently observed by the BaBar and Belle collaborations whose values 
are in good agreement with each other (Table I): the world average value is 

J+p = -0.109 ±0.020 . (27) 

The direct CP asymmetry data for the other B —>■ Ktt modes involve large uncertainties. 

We comment on the values of the ratios between the BRs for the B — > Kit modes, R±, 
R c , and R n , which are obtained from the experimental results given in Table I: 

Ri = = 0.82 ± 0.06 , (28) 

2B +0 

Rc = ^r = 1-00 ±0.09, (29) 
B + - 

Rn = jgoo = 0.79 ± 0.08 . (30) 



It has been claimed that within the SM, R c — R n « [1J, |3J|. From their definitions, it 
is indeed clear that R c ~ R n , if the r 2 -order terms including r EW or r c are negligible. In 
other words, any difference between R c and R n would arise from the contributions from the 
subdominant r 2 -order terms including r EW or r c . The above experimental data show the 
pattern R c > R n jlJ . 0], which would imply an enhancement of the electroweak penguin 
and/or the color-suppressed tree contributions. 
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TABLE I: Experimental data on the CP-averaged branching ratios {B in units of 10 -6 ), the direct 
CP asymmetries (Acp), and the mixing-induced CP asymmetry (S Ks7T o) for B — ► Kir modes. 
The Sk-zO is equal to sin(2<Aj ) in the case that tree amplitudes are neglected for B° — > K s ir° 



±iie (jf^^jrO is eui. 



9, 



mini 



CLEO 



Belle 



BaBar 



Average 



B(B ± - 




B{B ± - 




B{B° - 




B(B° - 





18.J 



12.9 



18.0 



12. 



+3.7+2.1 
-3.3-1.8 

-2.4+1.2 
-2.2-1.1 

-2.3+1.2 
-2.1-0.9 

-4.0+1.7 
-3.3-1.4 



22.0 ± 1.9 ± 1.1 
12.0 ± 1.3±J;| 

18.5 ±1.0 ±0.7 
11.7 ± 2.3±i;| 



26.0 ± 1.3 ± 1.0 
12.0 ±0.7 ±0.6 
17.9 ±0.9 ±0.7 
11.4 ±0.9 ±0.6 



24.1 ± 1.3 

12.1 ±0.8 

18.2 ±0.8 
11.5 ± 1.0 



A 0+ 
4+o 



A, 



A m 



CP 

00 
CP 



0.18 ±0.24 
-0.29 ±0.23 
-0.04 ±0.16 



0.05 ±0.05 ±0.01 
0.04 ±0.05 ±0.02 
-0.101 ±0.025 ±0.005 
-0.12 ±0.20 ±0.07 



-0.087 ±0.046 ±0.010 

0.06 ±0.06 ±0.01 
-0.133 ±0.030 ±0.009 
-0.06 ±0.18 ±0.06 



-0.020 ± 0.034 
0.04 ± 0.04 
-0.109 ±0.020 
-0.09 ±0.14 



S 



0.30 ±0.59 ±0.11 



0.35±§ |§ ± 0.04 



0.34 ±0.29 



We note that in this analysis we do not consider B — ► modes simultaneously with B — > 
Kit modes, though they can be connected to each other by using flavor SU(3) symmetry. 
The reason is that we do not want that our analysis to be spoiled by the unknown effects of 
the flavor SU(3) breaking relations after the inclusion of FSI effects. 

We remind that assuming the conventional hierarchy as in Eqs. (115)) and (|16j) . Aq° p is 
expected to be almost the same as Aq P : in particular, they would have the same sign. 
However, the data show that Aq P differs by 3.4<r from A^p. This is a very interesting 
observation with the new measurements of A^p by BaBar and Belle, even though the mea- 
surements of Agp still include sizable errors. One may need to explain on the theoretical 
basis how this feature can arise. 



B. Global x 2 Analysis and Theoretical Implications 

Based on the current experimental data shown in Table I, we critically investigate their 
implications to the underlying theory on the B —>■ processes. There are nine observables 
available for the B — > modes, as shown in Table QJ However, there are ten theoretical 
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parameters (\P\, r T , r EW , r c , 5 T , 5 EW , 5 C , cxi/2, «3/2 and 03 (= 7 of the unitary triangle)) 
relevant to the above nine observables, neglecting the very small terms of the annihilation 
contribution 1 r A , and the color suppressed electroweak penguin contribution . We have 
fixed sin20i(= (3) = 0.726 to its central value of the experimental measurements. Therefore, 
we have to fix at least one of the theoretical input parameters by assuming a model-calculated 
value or the previously measured central value in order to do x 2 analysis. 

1. Global x 2 Analysis within Isospin-phased Topological Approach 

TABLE II: x 2 fit for cases(a-d) and their combinations, including all sdSPs and IdSPs. Fixed 
parameters in each column appear in parentheses. 





(a) 


(b) 


(a+b) 


(c) 


(d) 


(c+d) 


Xmin/ d -°- f 


0.0/0 


0.0/0 


0.0/1 


0.1/0 


2.02/0 


2.04/1 


P 


22.7 ±2 .3 


23.56 ± o.78 


22.9-to.80 


21.9±i. 8 


24.0±i. 3 


23.8±i.3 


«l/2 


-0.03 ±0 .4i 


-0.10±o.25 


-0.05± .27 


0.52 ±0 .39 


0.25 ±0 .7i 


0.35 ± o.25 


"3/2 


0.23±o.7i 


0.0 ±0 .35 


0.18± .36 


-l.04io.68 


3.048±o.095 


3-04-to.087 


5t 


0.31±o.46 


0.32±o.i3 


0.28±o.n 


0.128±o.046 


0.20 ±0 .io 


0.19±o.n 


&EW 


l-68±o.69 


2.00 ±0 .34 


l-73±o.29 


(=*r) 


0.55 ±2 .o 


(=*r) 


Sc 


-3.01-to.43 


-2.90 ±0 .29 


-3.01± .35 


— 3.92-to.65 


(=5 P ) 


(=5 P ) 


r T 


0.19 ±0 .24 


(0.21) 


(0.21) 


0.58 ±0 .i8 


0.33±o.i4 


0.31 ±0 .ii 


r EW 


0.35 ±0 .i3 


0.322-to,080 


0.347±o.087 


0.01±o.33 


0.16 ± o.20 


0.142±o.04o 


r c 


0.28-to.24 


0.30 ±0 .i9 


0.29 ±0 .2i 


0.54 ±0 .28 


0.00 ±0 .i5 


0.00 ±0 .i5 


03(°) 


(60.0) 


50.9 ±5 .5 


(60.0) 


73.3-t8.4 


64.0±io.o 


64.8±io. 4 



We do the x 2 analysis to investigate sources of physics beyond the SM step-by-step as 
follows: for 9 observables (4 branching ratios, 4 direct CP asymmetries and one indirect CP 
asymmetry) with 9 input parameters by fixing case(a): 3 = 60° (but p oc /t|P| 2 |A t | 2 , oti/ 2 , 
a 3/2, &TEwc an d t tewc are ^ ree parameters.), or case(b): r T = 0.21 (the pQCD central 
value) or case(c): 5t = Sew or case(d): 5q = (i.e. 5' c = 5' P ). We further exercise by fixing 

1 After the annihilation term r A is neglected, the observable A^p still remains non-zero due to non-zero 
values of the isospin phases. In our approach, all nine observables still remain relevant. 
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case(e): the values of all sdSPs zero (typical isospin analysis), where we have three fewer 
parameters, or case(f): the values of all IdSPs zero (usual topological quark diagrammatic 
approach), in which we have two fewer parameters; and combinations of cases(a-f). 

In Table |H1 we show the results of x 2 -fhti n g of 9 observables by parameters including 
both IdSPs and sdSPs. For the minimum of x 2 , i n cases(a, b, c, a+b) Xmm vame are 
almost zero, whereas in cases(d) Xmm becomes larger than 2, which means that Sc = Sp 
for sdSPs is not a good assumption. From now on, we do not consider the case(d) any 
more, except the combination with case(c) which assumes Sc = Sp as well as St = Sew- 
We searched the whole range of parameter spaces for global Xmim an d we obtained the large 
IdSPs: a 1 / 2 in the range —0.9 < ai/ 2 < 0.5 (in the unit of radian) and «i/ 2 = is allowed 
within the range of la in cases(a, b, d, a+b); and for 0:3/2, in cases(a, b, a+b) 0:3/2 = 
is allowed within the la level; whereas in case(c) 0:3/2 = is not allowed in la. Central 
values of St are distributed in the range 0.7° < St < 18°, and central values of Sew are laxge 
and about 96° < Sew < 115° in cases(a, b, a+b); however, the assumption St = Sew is 
quite valid with larger r T , much larger than typical pQCD estimates. Central values of Sc 
are very large and almost 180° in all cases with the exception of cases(d, c+d) in which 
Sc — dp is assumed. For r T , in cases(a, b, a+b) central values are in good agreement 
with its pQCD value r T = 0.21. As for r EW and r c , in general large values are preferred for 
better \ 2 fitting; central values of r EW are around 0.35, much larger than its typical pQCD 
value r EW = 0.14, in cases(a, b, a+b). For r c , cases(a, b, c, a+b) exclude the pQCD 
value r c = 0.018 by more than la. 



TABLE III: CP asymmetries estimated from best-fit parameters for each case in TABLE UTI 



CP asym. (Exp.) 


(a) 


(b) 


(a+b) 


(c) 


(d) 


(c+d) 


A+p (-0.109 ±0 .o2o) 
A+P P (0.04±o.o 4 ) 
A° c + P (-0.020 ±0 .o34) 
A° C ° P (-0.09±o.u) 

Sk s tt° (0-34±0.29) 


-0.109±o.28 
0.04 ±0 .34 

-0.020 ± o.h 

-0.09±o.28 

0.34 ±0 .38 


-0.109 ±0 .067 

0.04 ±0 .i7 

-0.020±o.055 

-0.09±o.i6 

0.34 ±0 .26 


-0.109±o.o80 

0.04 ±0 .2i 
-0.020±o.o7o 
— 0.08±o.2o 

0-33±o.32 


-0.107±o.33 
0.04 ±0 .48 

-0.021 ±0 .30 

~0.09±o.57 

0.43 ±0 .30 


-0.106 ±0 .ioo 
0.026 ±0 .o85 

-0.014 ±0 .122 

-0.18±o.2i 

0.09 ±0 .28 


-0.106 ±0 .038 

0.025 ±0 .o43 

— 0.013±o.O53 

-0.19 ± o.n 
0-12±o.o5 



In Table IILT1 we show the estimated CP asymmeries from the best fit parameters for each 
case in Table |Hj to be compared with the experimental average values of Table [IJ We do 



12 



not show the estimated branching ratios because there are almost no differences between 
the experimental values and the theoretically estimated values from the best fit p arameters. 

A few comments are in order here: (i) Contrary to the previous finding jlj . within 
the topological approach, case(c) confirms 5t = Sew of Ref. j^. It can be understood by 
the fact that the topology of the tree diagram is quite similar to that of the EW penguin 
diagram, (ii) Case(c+d) is very interesting since it gives the pQCD estimated values on 
the parameters, r T , r EW , r c and 03 (see Eq. (HHJ)). That is, if we include FSI, we can 
fit the parameters within the SM without invoking any unknown new physics effects at 
all. (iii) The IdSPs from FSI can be quite large, and therefore, we cannot ignore the final 
state re-scattering effects even in B meson decays, (iv) Indeed, it has been long advocated 
that charming-penguin contributions can increase significantly the B — > Kit rates and yield 
better agreement with experiment 13211 . Strong phase S c from charming-penguin has been also 
bund substantially large, 5 C ~ 20° [33]. (v) Long distance strong phases has been modelled 
2^ | as re-scattering of some intermediate two-body states with one particle exchange in the 
t channel and the absorptive part of the re-scattering amplitude via optical theorem. Large 
long distance contributions are found as 5t — 8c ~ 90° with FSI, too. (vi) It has been 
found 3^ that by using general features of soft strong interactions soft scattering does not 
decrease for large mg, and inelastic processes are expected to be leading sources of strong 
phases 2 . Please note that in this paper we only include the elastic FSI contributions through 
isospin phases. 



2. Null Hypothesis Test by Statistical p— value 

Here we make a hypothetical test by using a statistical p— value 3 , which is defined by the 
integral J°2 g(t;d)dt, where Xmin * s the minimum \ 2 value obtained and g(t;d) the proba- 

'Mil in 

bility density function for x 2 distribution with given d degree of freedom. One hypothesis 

2 ^From the viewpoint of Regge theory, it is pointed out 34] that even in the heavy quark limit there are 
non-vanishing inelastic non-perturbative FSI effect remaining on each individual channels. However, in 
[s^l authors showed the inelastic FSI effects will be cancelled by each other when all individual channels 
are summed up. Because in B hadronic decays the 6-quark mass is finite, therefore, such a cancellation 
might be incomplete and the non-perturbative inelastic effects could be sizable ^tf^. 

3 In the sense of statistics, the p— value is the probability that the observed value happens to be observed 
under the "null hypothesis" . When the p— value is smaller than a threshold (usually 5% is used) , it can 
be said that the null hypothesis is disfavored and the "alternative hypothesis" is favored. 
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TABLE IV: The x 2 -result with fixed r's (i.e. r T = 0.21, r EW = 0.14 and r c = 0.018). In "free" 
case, p, a.i/2, «3/2) $t,ew,C and (ps are all free parameters. In each cases 2 IdSPs, 4 branching 
ratio, 4 direct and 1 indirect CP assymmetries at the x 2 -minimum are also shown. 



Strong Phases 


both IdSPs and sdSPs 


only 


sdSPs 


only IdSPs 




(a) + free </> 3 


(c+d) 


(f) 


(f+c) 


(e) 


Xmin/ d -°-f 


2.3/2 


4.3/3 


8.4/4 


12.8/5 


16.2/5 


p— value(%) 


32 


23 


7.8 


2.5 


0.6 



which we want to test is the assumption that the parameters ry, r EW and rc take their 
pQCD estimated values, e.g. r T ~ 0.21, r EW ~ 0.14 and r c ~ 0.018. We call this the 
"amplitude assumption". The other hypothesis is the relation among sdSPs: 5t = Sew & n d 
5c = Sp, which we call the "strong phase assumption". 

In Table IIV1 we summarize the Xmm an d its p— value with amplitude and strong phase 
assumptions for each case — case with only IdSPs, case with only sdSPs and case with 
both IdSPs and sdSPs. In the case with both IdSPs and sdSPs, p— value is about 32% 
with non-constrained sdSPs (free), and one with constrained sdSPs (case(c+d)) is still 
around 23%. In the case with only sdSPs (cases(f, f+c)), p— value in (f) is 7.8%, while 
in (f+c) the probability is reduced to 2.5% and this case is almost ruled out. In the case 
with only IdSPs (case(e)), the p— value is only 0.6% and this case must be completely ruled 
out. Therefore, this test confirms our previous claim that the generalized "isospin-phased 
topological approach" can fit the B — > decays well with the SM values of the parameters, 
which are estimated from pQCD, without invoking any unknown new physics effects. We 
again conclude that we cannot ignore the final state re-scattering effects and IdSPs even in 
B meson decays. 

3. Dependence of Sk s -k and A^p on r EW and r c 

Now we would like to make a few comments on the sensitivity of the observable Sk S w° to 
the parameter r c . The experimental value of S^gn is much smaller than that of Sj/^k s - ^ s 
implied by (J26|) . the theoretical prediction of Sk S tt° can be very sensitive to the parameter 
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TABLE V: \ 2 result for (A) varying S Ksn o = 0.2, 0.4, 0.6, 0.7, and, (B) varying A+° p = 
+0.04, +0.02, —0.02, —0.04. Error is assumed to be 20%. Concerning the input parameters, we 
choose case(a) (i.e. IdSPs and sdSPs are included, fa = 60° is fixed, and sin 2/3 = 0.726 is used). 



(A) 



S Ks7T o 


0.20+0.04 


0.40+0.08 


0.60+0.12 


0.70+0.14 


x 2 - 

Amm 


0.0 


0.0 


0.0 


0.0 


r 

' EW 


0.34+0.16 


0.36+0.12 


0.36+0.10 


0.36+0.10 


r c 


0.39+0.09 


0.23+0.09 


0.08+0.07 


0.06+0.11 



(B) 





+0.04 ± 0.008 


+0.02 ± 0.004 


-0.02 + 0.004 


-0.04 + 0.008 


x 2 - 


0.0 


0.0 


0.0 


0.0 


T 

' EW 


0.35 + 0.13 


0.30 + 0.16 


0.38 + 0.16 


0.38 + 0.10 


r c 


0.28 + 0.25 


0.28 + 0.24 


0.15 + 0.09 


0.26 + 0.26 



r c . For illustration, in part (A) of Table IVl we vary Sk s7 t° around the present experimental 
value, specifically, S Ksn o is assumed to be (0.20 ± 0.04) , (0.40 ± 0.08) , (0.60 ± 0.12) and 
(0.70 ± 0.14), respectively. Here just for the illustrative purpose, we set 20% errors in each 
case. (Also, to be consistent, we set 20% errors to all the data whose current errors are 
larger than 20%, such as A^p.) It can be seen that as Sk s it° varies, the allowed region for 
r c varies in a clearly correlated manner. Just for comparison, in (B) of Table 13 we also 
present the case when the value of Aq P varies. We recall that the averaged experimental 
value of Agp is positive, in contrast to the value of A^p. Again for the illustrative purpose, 
A£p is assumed to be (±0.04 ± 0.008), (±0.02 ± 0.004), respectively. (To be consistent, we 
also set 20% errors to all the data whose current errors are larger than 20%, such as Sk S tt° 
and Aq° p .) In contrast to case (A), the dependence of r EW and r c on A^ p is very obscure. 

IV. CONCLUSIONS 

All CP-averaged branching ratios of B — > Kir modes have been recently measured within 
a few percent errors by the BaBar and Belle collaborations. The observations of the direct 

15 



CP asymmetry in B° — > K ± tt t have also been recently achieved at the 5.7a level by BaBar 
and Belle. Certain experimental data (e.g., the branching ratios (BRs)) for B — > Kn are 
currently more precise than the theoretical model predictions based on QCD factorization 
or perturbative QCD. Thus, with the currently available precision data, it is very important 
to investigate these modes as generally and critically as possible. The importance of the 
soft final state re-scattering effects, especially, in a model-independent topological quark 
diagrammatic analysis of B — > Dit, B — > Kit and B — > txtt decays, has also been recognized 
recently. It is certainly conceivable that final state interactions can modify the predictions 
based on the short distance diagrammatic analysis. In this work, we studied the decay 
processes B — > Kir by using generalized decay amplitudes within the so-called "isospin- 
phased topological approach" , i.e. with both the strong phases of the final state re-scattering 
and the strong phases of the topological quark diagrammatic origin. We are mainly interested 
in investigating whether the conventional SM predictions are consistent with the current data 
even after we include final state interaction effects. In order to find the most likely values 
of the parameters in a statically reliable way, we used the \ 2 analysis. 

Our result shows that: (i) Contrary to the previous finding within the topological 
approach, we confirmed 5t = Sew of Ref. which can be understood by the fact that 
the topology of the tree diagram is quite similar to that of the EW penguin diagram, (ii) 
If we include final state interaction effects, we can fit theoretical input parameters within 
the SM without invoking any unknown new physics effects at all. (iii) The final state 
re-scattering phases can be quite large, and therefore, we cannot ignore the final state 
re-scattering effects even in B meson decays, (iv) We also found that there are strong 
correlations between the parameter r c and the time dependent indirect CP observable 
S Kg7T o: If S Ks7T o « Sj/^k s , then r c < 0.1. However, the present experimental value, 
Sk s -k° <S Sj/^k 3 , directly implies very large r c like 0.3 ~ 0.4, if final state re-scattering 
phases are not considered. 
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